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Qh Abstract 

We study a supersymmetry breaking mechanism in the context of a minimal 
anomalous extension of the MSSM. The anomaly cancellation mechanism is achieved 
through suitable counterterms in the effective action, i.e. Green-Schwarz terms. 
We assume that the standard MSSM superpotential is perturbatively realized, i.e. 
I/-) all terms allowed by gauge symmetries, except for the /it-term which has a non- 

perturbative origin. The presence of this term is expected in many intersecting D- 
brane models which can be considered as the ultraviolet completion of our model. 
We show how soft supersymmetry breaking terms arise in this framework and we 
study the effect of some phenomenological constraints on this scenario. 
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1 Introduction 



The LHC era has begun and the high energy physics community is analyzing and dis- 
cussing the first results. One of the key goals of LHC, besides shedding light on the 
electroweak (EW) symmetry breaking sector of the Standard Model (SM), is to find some 
signature of physics beyond the SM. Supersymmetric particles and extra neutral gauge 
bosons Z 1 are widely studied examples of such signatures. A large class of phenomenologi- 
cal and string models aiming to describe the low energy physics accessible to LHC predict 
the existence of additional abelian U(l) gauge groups as well as iV = 1 supersymmetry 
softly broken roughly at the TeV scale. In particular in string theory the presence of 
extra anomalous £/(l)'s seems ubiquitous. D-brane models in orientifold vacua contain 
several abelian factors and they are typically anomalous [1]-[16]. In [17] we studied a 
string inspired extension of the (Minimal Supersymmetric SM) MSSM with an additional 
anomalous U(l) (see [18] for other anomalous U(l) extensions of the SM and see [19] 
for extensions of the MSSM). The term anomalous refers to the peculiar mechanism of 
gauge anomaly cancellation [20] which does not rely on the fermion charges but rather 
on the presence of suitable counterterms in the effective action. These terms are usually 
dubbed as Green-Schwarz (GS) [18, 21] and Generalized Chern-Simons (GCS) [22]-[26]. 
They can be considered as the low energy remnants of the higher dimensional anomaly 
cancellation mechanism in string theory. In our model we assumed the usual MSSM su- 
perpotential and soft supersymmetry breaking terms allowed by the symmetries (the well 
known result [27]). In this paper we address the question of the origin of the latter in the 
context of a global supersymmetry breaking mechanism. This means that we do not rely 
on a supergravity origin of the soft terms but rather on a local setup based for example 
on intersecting D-brane constructions in superstring theory in which gravity is essentially 
decoupled (see for instance [28] for a recent attempt in this direction). Moreover in [17] 
we made the assumption that all the MSSM superpotential terms were perturbatively 
realized, i.e. allowed by the extra abelian U(l) symmetries. In the following we assume 
instead that the /x-term is perturbatively forbidden. The origin of this term is rather non- 
perturbative and can be associated to an exotic instanton contribution which naturally 
arises from euclidean D-brane in the framework of a type IIA intersecting brane model 
(see [29] and references therein). 

The paper is organized as follows: in Sec. 2 we describe the basic setup of the model and 
we discuss the perturbative and non-perturbative origin of the superpotential terms. We 
argue how the latter can naturally come from an intersecting D-brane model considered 
as the ultraviolet (UV) completion of our model. In Sec. 3 we describe the (global) 
supersymmetry breaking mechanism that gives mass to all the soft terms. In Sec. 5 we 
compute the gauge vector boson masses while in Sec. 4 we study the scalar potential of 
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the theory in the neutral sector. In Sec. 6 we describe the neutralino sector while in 
Sec. 7 we describe the sfermion mass matrices. In Sec. 8 we study the phenomenology of 
our model and the bounds that can be put by some experimental constraints. Finally in 
Sec. 9 we draw our conclusions. 

2 Model Setup 

The model is an extension of the MSSM with two extra abelian gauge groups, U(1)a and 
U(1)b- The first one is anomalous while the second one is anomaly free. This assumption 
is quite generic since in models with several anomalous U(l) symmetries there exists a 
unique linear combination which is anomalous while the other combinations are anomaly 
free. The charge assignment for the chiral superfields is shown in Table 1. The vector 
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Table 1: Charge assignment. 



and matter chiral multiplets undergo the usual gauge transformations 

V -> V + i(A-A jf ) 

$ ->■ e- iqA <5> (1) 

The anomaly cancellation of the U(1)a gauge group is achieved by the four dimensional 
analogue of the higher dimensional GS mechanism which involves the Stiickelberg super- 
field S = s + Idips + 9 2 Fs transforming as a shift 

S ->■ S- 2iMy A A (2) 

where M Va is a mass parameter related to the anomalous U(1)a gauge boson mass. It 
turns out that not all the anomalies can be cancelled in this way. In particular the so called 
mixed anomalies between anomalous and non anomalous £/(l)'s require the presence of 
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trilinear GCS counterterms. For further details about the anomaly cancellation mecha- 
nism see Appendix A (see also for instance [17] and [24]). The effective superpotential of 
our model at the scale E = M Va is given by 

W = W M ssm + \e~ kS H u H d + m$+$- (3) 

where W M ssm is given by 

W M ssm = yHQiU?H u - y l3 d Q % D]H d - y^UE]H d (4) 

which is the usual MSSM superpotential without the //-term which is forbidden for a 
generic choice of the charges qn u and qn d . The second term in (3) is the only gauge 
invariant coupling allowed between the Stiickelberg superfield and the two Higgs fields. 
This is the only allowed coupling with matter fields for a field transforming as (2). We will 
argue later about how non perturbative effects can generate such a term. The last term 
in (3) is a mass term for $ ± which are charged under both U(1)a and C/(1) B . These fields 
have been considered as supersymmetry breaking mediators in the context of anomalous 
models by Dvali and Pomarol [30]. They play a key role in generating gaugino masses. 
In the effective lagrangian, besides the usual kinetic terms (they are charged under both 
U(1)a and U(1) B ), the two U(l) B fields $ ± couple to the gauge field strength W° through 
the dimension six effective operator 

$+$- 

Cg = Ca-gj-WZWaa (5) 

where a = A, B, Y, 2, 3, A is the cut-off scale of the theory while c a are constants that 
have to be computed in the UV completion of the theory. 

The non perturbative term in (3) is expected to be generated in the effective action 
of intersecting D-brane models which can be considered as the UV completion of our 
model. This is the leading order term when the coupling H u H d is not allowed by gauge 
invariance. In string theory there are many axions related to the GS mechanism of 
anomaly cancellation which are charged under some Ramond-Ramond (RR) form. For 
example in type IIA orientifold model with D6-branes, axion fields are associated to the 
C 3 RR-form (see for a recent review [31]). Instantons charged under this RR-form, such as 
euclidean E2-branes wrapping some 73 3-cycle in the Calabi-Yau (CY) compactification 
manifold, give a contribution to the holomorphic couplings in the N — 1 superpotential. 
Our analysis does not rely on any concrete intersecting brane model but rather on the 
generic appearance of such instanton induced terms. The exponential suppression factor 
of the classical instanton action is 

e -Voi E 2/ss (6) 

where Vo1e2 is the volume of the 3-cycle in the CY wrapped by a i?2-brane measured in 
string units while g s is the string coupling. Such exponential factor is independent from 
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the d — 4 gauge coupling and thus this instanton is usually termed as stringy or exotic 
instanton (see [29] and [32] and references therein). Moreover the instanton contribution 
can be sizable even in the case g s <C 1 if Vo1e2 <C 1 measured in string units. 

In type IIA orientifold models with intersecting branes the complexified moduli, whose 
imaginary part are the generalized axion fields (depending on the cycle 73), can be written 

as 



Ui = e- V / n 3 + t / C 3 (7) 

where <p is the dilaton, Q3 is the CY volume 3-form (which is a complex form) and C3 
is the RR-form. The integral of this form is dual to the axion whose shift symmetry is 
gauged in the GS mechanism. The generic contribution of an E2 instanton is formally 
given by 

n 

1=1 

where $ ai ,&; are chiral superfields localized at the intersection of two D6-branes described 
by open strings while Se2 denotes the instanton classical action 



e Se2 = exp 



2tt / 1 
% \9sj- 



Re(Q 3 ) -i Q 



(9) 



This result can be immediately extended to the supersymmetric case which involves the 
complete Stiickelberg multiplet. The appearance of the exponential suppression factor 
is dictated by the fact that the superpotential is a holomorphic quantity. Thus the only 
allowed functional dependence on the string coupling g s = e <<f> and the axionic field is an 
exponential. Any other dependence can be excluded due to the shift transformation (2). 



3 Supersymmetry Breaking 

The D-term contribution of the U(1)a vector multiplet Va relevant to supersymmetry 
breaking is given, in the limit of vanishing kinetic mixing 8ya, $ab = 0, by the following 
lagrangian: 

C = l -D A D A + 9Aqi<f>\ D A<t>i + iD A (10) 

i 

where the sum is extended to all the scalars charged under the U(1)a- There is no D- 
term contribution related to the U(1)b except that of ± since all the MSSM chiral fields 
are uncharged under U(1)b (see Table 1). The last term in (10) is a tree-level field de- 
pendent Fayet-Iliopoulos (FI) term which comes from the supersymmetrized Stiickelberg 
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lagrangian 

1 



t (S + & + 2M VA V A y 



+ 

e 2 e 2 



= M Va (S + tf)v A \ e2§2 + . . . 

= M VA aD A + ... (11) 

where in the last line a denotes the real part of the lowest component of the Stiickelberg 
chiral multiplet s = a + itp. The fields a and ip are called the saxion and the axion 
respectively 3 . We assume that the real part a gets an expectation value. This gives 
a contribution to the gauge coupling constants which can be absorbed in the following 
redefinition 

! ! l b aa (a) (12) 



16#2 Ta lQg 2 a r a 2" 

where the gauge factors r a take the values 1,1,1,1/2, 1/2 and the b aa constants are given 
in (101). The tree-level FI term is then given by 

i = M Va (a) (13) 

Moreover in the following we assume that 1-loop FI terms are absent (see the discussion 
in [33]). The FI term induces a mass term for the scalars. This can be seen by solving 
the equations of motion for D A 

D A + J29Aqi4<f>i + Z = U (14) 

i 

where the index i runs over all chiral superfields. The D-term contribution to the scalar 
potential is given by 



l^ + gA^M^j ( 15 ) 



The quadratic part gives the scalar mass term 

^2^9Aqi\<Pi\ 2 = ^2m 2 i \<f) i \ 2 (16) 



where we have defined 



with 



mf = £g A qi = (a) g A M VA qi = qiin 2 (17) 



m\ = {a)g A M VA =g A i (18) 

The typical scale for the mass is of the order of few hundreds of GeV if M Va ~ 
(a) ~ 1 TeV and g A ~ 0.1. It is interesting to note that in this scenario a low subTeV 



3 with a slight abuse of notation with respect to the previous section where we denoted the dilaton 
with (p. 
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supersymmetry breaking scale is due to the Stiickelberg mechanism which gives mass 
to Va- This is the most important difference with the scenario proposed in [30], where 
the scale is dynamically generated by some dynamics in a strong coupling regime. 

Mass terms for the gauginos, i.e. A a A a , are generated by the dimension six effective 
operator (5) in the broken phase where ± get vacuum expectation value (vev). The 
contribution coming from this mechanism is 

M _ c + _ m(v% + vj) 

lvla ~ Ca ~ Ca 2 A 2 

where v±/ \/2 = (</>±) and where in the right hand side we have used the F-term equations 
of motion for F ± 

dW* 

F ± = -^i- = -m^* (20) 

having assumed m real without any loss in generality. We assume c a = c for each a. This 
is an assumption of universality as a boundary condition at the cutoff scale A which does 
not affect in a crucial way our analysis. In section 4 we study the scalar potential of our 
model and we derive the conditions for having a vev for ± different from zero. Since we 
are breaking supersymmetry in the global limit in which the Planck mass Mp — > oo the 
F-term induced contribution to the scalar masses 



™2 ( F ± ,, rn 

M 2 



p 

vanishes leaving (17) as the leading contribution. 

The requirement of gauge invariance of the superpotential implies the following con- 
straints on the U(1)a charges 

Que = -qq - Qh u 
Qd- = -qq - qu d 

qsc = -qL-qH d (22) 

and 

k = qHu + qHd (23) 
2M Va v ; 

As we said at the beginning of this section we assume that the net kinetic mixing between 

U(1)y and U(1)a vanishes 4 . There are two contributions for the U(l)y — U(1)a kinetic 

mixing: the 1-loop mixing 5 Y a and b YA coming from the GS coupling SW y Wa (see 

eq. (100)). The following conditions imply a bound on the charges 

5 YA = ^q f Yf = 
f 

b YA = ^ g Jr / = (24) 

/ 



4 We postpone the discussion about the kinetic mixing between U(1) A and U(1)b to the next section. 
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where the sum is extended over all the chiral fermions in the theory. The con- 
straints (24) can be solved in terms of q Q and q L . By using the conditions (22) we 
get 

1l = \ (3<1h u ~ 4:q Hd ) 

QQ = -^(^H u -2q Hd ) (25) 

The positive squared mass condition for the sfermions 

m) = g A q f M VA (a) > (26) 

implies qf > for all the sfermions having assumed without loss of generality (a) > 0. 
Using the constraints (22) and (25) we get the allowed parameter space 

5 3 
Qh u < , -q Hu < q Hi < ^q Hu (27) 



4 Scalar Potential 



The key ingredient in our model is the instanton induced term in (3) which couples the 
Stiickelberg field to the Higgs fields. The 9 2 component of this superpotential term gives 
the following contribution to the lagrangian 

Wi ns t\e2 = Xe~ kS H u H d \ d 2 

= Xe~ ks h u F d + Xe- ks F u h d - Xke~ ks F s h u h d + 

V2Xe~ ks k (h u i> s h d + M^«) - Xe- ks k 2 h u h d i> s i>s (28) 

where F Uyd are the F-terms of H U)d . Solving the F-terms equations for H u and H d we get 
the following contributions for the instanton induced term in the scalar potential 

V mst = 2X 2 e- 2ka hlK + 2X 2 e~ 2ka h\h d + Xke- ka (e~^F s h u h d + h.c.) (29) 

In the following we assume that a gets a vev different from zero and that the mass of 
this field is much higher than A so that its dynamics is not described by the low energy 
effective action. From the point of view of the UV completion (for example a type IIA 
intersecting brane model) this amounts to saying that the closed string modulus related to 
a is stabilized. Moreover we made the assumption that the same dynamics that stabilizes 
a also fixes F s . By supersymmetry the saxion field a, being part of the Stuckelberg 
multiplet, has a tree-level mass My A . Thus if we want to consider a frozen dynamics 
for a at the TeV scale we have to assume a mass parameter for the anomalous U{1)a 
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just slightly above the TeV scale, i.e. M Va > 1 TeV. In this way the effective instanton 
induced potential at a scale E ~ 1 TeV is thus given by 

V mst = 2\ 2 e- 2k ^hlh u + 2\ 2 e~ 2k ^h\h d + \ke~ k ^ ((F s ) e -^h u h d + h.c.) (30) 

The first two terms are /i-terms while the third one is a b-term. The complete effective 
scalar potential is given by 

v = (H 2 + 0(I^T + l^r) + (H 2 + 0(l^l 2 + l^l 2 ) 

+ (|m| 2 + mj + )|0+| 2 + (H 2 + ml-)\<p-\ 2 ) 
+ [be-** {Kh~ d - h°X) + h.c] 

+ \{9l + 9 2 v) (K\ 2 + \K\ 2 - K\ 2 - \h-\ 2 ) 2 + \g\ \h#% + 



"d 



+\g\ (K\ 2 + \K\ 2 ) + qHd {\h d \ 2 + \h d \ 2 ) + |0+| 2 - |0-| 2 ] 2 

+^io^ + i 2 -n 2 ] 2 (si) 



where 



At = v / 2Ae- fe<Q) (32) 
6 = \ke- k ^(F s ) (33) 

These relations give a solution of the well known /i-problem since both terms have a 
common origin (see the analysis in Sec. (8.2)). The soft squared masses are generated by 
the FI U(1) A term 

™Z U = <lHjn\ (34) 

m l d = QH d ml (35) 

m l+ = m \ (36) 

m\- = -m\ (37) 

with m| given by (17). The scalar potential depends on the following new parameters: 
(a), (F s ), A, m, g AB , qn uA , %• 

In order to have a vacuum preserving the electromagnetism the charged field vevs 
must vanish. Thus we are left with the problem of finding a minimum for the neutral 
scalar potential 

v = (H 2 +<)i^r+(i/ii 2 +<)i^r-(^^«+^.) (38) 

+ (|m| 2 + mj + )|0 + | 2 + (H 2 + ml-)\(f)'\ 2 
+ \{9l + 9 2 y){\hl?-K\ 2 ) 2 

+\g\ (^j^i 2 +^i^r+i0 + i 2 -n 2 ) 2 



Since there are no D-flat directions along which the quartic part vanishes, the potential is 
always bounded from below. To find the minimum we solve dV^/dz 1 = where the scalar 
field z % runs over {ip, h° u , h Q d , </> + , </>"}. The conditions for having a non-trivial minimum 
boils down to the same condition of the MSSM 

& 2 >(H 2 + 0(H 2 + <) (39) 

Moreover in order to generate a mass term for the gauginos (see eq. (19)) the condition 
v - ^ must hold since v + = due to the positive sign of the coefficient of the <p + 
quadratic term in (36). This implies the following condition for the coefficient of the <p~ 
quadratic term 

H 2 + mj_ < (40) 

The minimum is attained at ip = <p + = 0. Actually since the potential for the axion (p 
is periodic the minimum condition holds for ip = 2nir/k with n G Z. All these minima 
are physically equivalent and thus we arbitrarily choose n = 0. The remaining three 
conditions imply the following constraints on the parameters 

m l d + ^ -btp + -(gy + gl)v 2 c 2 p + ^g 2 A q Hd [v 2 (q Hd 4 + q Hu sl) - v 2 _] = (41) 

m\ u +fI 2 -bt/- l -{g 2 Y + gl)v 2 c 2p + l -g 2 A q Hu [v 2 {q Hd 4 + q Hu s 2 p) ~ v 2 _] = (42) 

(9a + 9%) vi - g\v 2 (q Hd c 2 + q Hu s 2 ) + 2 (\m\ 2 + mj_) = (43) 

where we have defined in order to keep a compact notation 

op = cos/3, S/3 = sm(3, tp = tan/3, c 2 p = cos(2/3), s 2 p = sin(2/3) (44) 

and as usual as tan/3 = v u /v<i- 

In the previous discussion we treated the scalar potential in an exact way. In the 
following we want to introduce some useful approximation in order to compute the mass 
eigenstates. Let us go back to the minima equations (41-43). Supposing 

v < v_ (45) 

we can neglect all the g A v terms. With this approximation the minima equations read 

m 2 hd + fi 2 -btp + ^(g 2 Y + g 2 2 )v 2 c 2 p = (46) 

m 2 K + fi 2 -bt- 1 - l -tf Y + g\) v 2 c 2 p = (47) 
{g 2 A + g 2 B )vl+2{\m\ 2 + ml-) = (48) 



where we have defined 



m h d = m h d ~ ^9 A qu d v^ (49) 
1 

T 



2' 

™L = m lu - n9 2 A qH u V 2 - (50) 
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Equations (46) and (47) have the same functional form as in the MSSM case. Moreover 
V- does not depend on any parameter of the visible sector. Within this approximation 
the dynamics of the fields ± is decoupled from that of the Higgs sector and thus the 
Higgs potential can be studied by fixing ± at their vevs. We get 

v - (H 2 +o(i^i 2 +i^r) + (H 2 +o(i^r+i^-r) 

+ [be~ ik * (h+h d - h°X) + h.c] 

+\{gl + a 2 y) (K\ 2 + K\ 2 - \h°f - |/^| 2 ) 2 + \ 9 l \hM* + hlh?\ 2 



1 2 



<iH u (K\ 2 + \K\ 2 ) + qHd (\h d \ 2 + \h- d \ 2 )-l^_ 



(51) 



neglecting further constant terms in t>_. Close to the minima the relevant term in the 
last line of eq. (51) is the double product of the Higgs part with the v 2 _ term. Hence by 
using (45) we finally get 

~ (H a + mL)(Kl 2 + l^l 2 ) + (H 2 + <)(W + l^| 2 ) (52) 
+ [be~ ik * (h + u h- d - h°X) + h.c] 

+ \(9 2 2 + 9 2 y) + \K\ 2 - W - |^-| 2 ) 2 + \gl \Kh d * + hlh~ d *\ 2 

This potential has the same form (except for the contribution of the exponential term in 
(p) of the MSSM potential and the corresponding minima equations are exactly given in 
eqs (46) and (47). Thus all the well known MSSM results apply here [36]. 

In particular one of the constraints is tp > 1.2 [36] which implies 5 m\ u < rh\ d - By 
using the equations (49) and (50) we get 

9aQh u (^M VA (a) - ^9av 2 ^J < g A qH d (^M VA (a) - ^av 2 ^ (53) 

By assuming M Va > 1 TeV, t>_ in the TeV range, qa ~ O(0. 1) the term between brackets 
is positive and we get the following constraint 

Qh u < q Hd (54) 

for the U(1)a Higgs charges. 



4.1 Higgs mass matrices 

We discuss the mass eigenvalues starting from the exact form of the scalar potential (31), 
switching to the approximated expression (52) when needed. In the neutral sector the 

5 The presence of the extra field ip does not affect this result since the minima conditions are the same 
as the MSSM. 
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singlet scalar + does not mix with any other scalar so it is a mass eigenstate with square 
mass 

M| + = 2\m\ 2 (55) 

The same holds for the imaginary part of (f>~ which becomes the longitudinal mode of the 
gauge vector Z 2 . The mass matrix for the real scalar fields {</?, im(/i°), Im(h^)} is given 
by 

/ btp ... ... \ 

M ( t l) = b bt^ 1 ... (56) 

y —bkv S/3 —bkvcp bk 2 v 2 cpsp J 

The determinant of this matrix is zero. Two eigenvalues are zero which correspond to 
the Goldstone modes of Z and Z\. The physical massive state is an axi-higgs state with 
mass given by 

2b 



S2/3 



1 ((1h u + q Hd ) 2 v 2 2 



1 



16 



'2/3 



(57) 



where we used the relation (23). The mass matrix for the real scalar fields {i?e(/i°), 
Re(h®), 0^ = Re(4>~)} reads as 



Mf e) 



(58) 



( (b 2 EW+ 9 A (lH d 2 )v 2 c} + btp 

~ b ~ [\9ew ~ 9 2 A qH d qH u ) v 2 c fi Sf i (\g 2 EW + g\q Hu 2 ) v 2 s} + bt' 1 



V 



-9 A qH d vv-Cp 



-9aQh u vvsi3 



(9 2 a + 9 2 b)v 2 



where g 2 EW = {g\ + g 2 ). The matrix can be diagonalized exactly but the results are cum- 
bersome and difficult to read. It is much more convenient starting from the approximated 
potential (52) neglecting the mixing between Higgses and <p~ . In this case we can apply 
the MSSM equations and get the following mass eigenvalues 



M, 



h°,H° 



Mi- ~ 

<VR 



1 

2 



2b 

S2f3 




(9 2 y + 9l> 2 + 2b(g 2 Y + g 2 )v 2 s w 



{9 2 a + 9 2 b) 



(59) 
(60) 



The charged sector is unchanged with respect to the MSSM, so 

2b 

S2/3 



M 2 H± 



+ M 2 ? 



(61) 



As in the standard MSSM case the mass of the lightest Higgs M h o has a theoretical 
bound. It is a well known problem in the MSSM that the upper bound [37] is not 
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compatible with the LEP bound [38]. In our case the bound is increased due to the 
presence of D^-term corrections 

Ml, < \{g 2 Y + gl)v 2 c% + ^g 2 A v 2 [q Hd + q Hu + {q Hd - q Hu ) c 2/3 f (62) 

where the first term is the MSSM bound. In principle, for arbitrary high values of gAqH d , 
gAqH u we get an increasing upper bound. However, as in the standard MSSM case, M% 
undergoes to relatively drastic quantum corrections [36]. Hence in Section 8 we consider 
tree-level masses for all the particles except for h for which we use the 1-loop corrected 
expression (see eq. (99)). 



5 Vector mass matrix 

We now discuss the vector mass matrix. All the neutral scalars could in principle take a 
vev different from zero, hence we assume 



v± 



( h °u,d) = 

The neutral vector square mass matrix in the base (Vb, Va, Vy, V 2 3 ) is 



M 



v 



( g\v\ 

9a9bvI g\ [{4q Hd 2 + s}q Hu 2 ) v 2 + v 2 ] + M, 

o \gAgYqH(P)v 2 
\ o 



2 

V A 



-\gAg2qH{fi)v z 



where 



(63) 
(64) 



(65) 



-\gYg2v 2 \g\v 2 j 



v A = 



V 



2 1 2 
v + + v_ 

2 1 2 

v u + v d 



(66) 
(67) 
(68) 



By taking M Va > 1 TeV (see Sec. 4), Va can be considered as decoupled from the low 
energy gauge sector (namely E < 1 TeV), and we can ignore with very good approximation 
any mixing term 6 involving Va- From now on we will apply this approximation. 



6 The kinetic mixing between ^(l)^ and U(1)b deserves some comment, in particular if we relax 
the My A > 1 TeV assumption. Actually the presence of this mixing turns out to be irrelevant for the 
phenomenology of the visible sector. Anyway one has to take into account that for Tr (qaQb) 7^ such a 
mixing arises at the 1-loop level. In such a case it can be assumed that the two 17(1) 's are in the kinetic 
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Since Vb is a hidden gauge boson, it is decoupled from the SM sector. The charged 
vector sector is unchanged with respect to the MSSM, so 



W ± = 



1 2 2 



(69) 
(70) 



6 Neutralinos 



In comparison with the standard MSSM we now have five new neutral fermionic fields: 
■ips, Aa, As, ± . However under the assumption My A > 1 TeV, ips and A a are not in the 
low energy sector because of the M Va mass term 7 . Thus we have 



Aieutralino mass ~~ o (^°) T ^N t l jV + ll - c - 



where 



(^°) T = (A S , 0", 0+, Ay, A°, hi hi) 
In this basis the neutralino mass matrix M.^ is written as 



(71) 



(72) 



/ 



N 



\ 



M B 
-g B V- 















-m 















Mi 


givd 

2 

givu 

2 



M 2 

92V d 
2 

92Vu 







(73) 



/ 



where \x is given in eq. (32). We remind that gaugino masses arise from the Dvali-Pomarol 
term (5). 

Ai factorizes in a 4 x 4 MSSM block in the lower right corner, and in a 3 x 3 new 
sector block in the upper left corner. The new sector block is given by the A# and ^ 
contributions. This last block has a MSSM-like structure that can be easily understood 



diagonalized basis with Tr (qaQb) = thanks to some additional heavy chiral multiplet charged under 
both U(1)a and U(1)b- These multiplets generate a counterterm in the effective theory that cancels 
against Sab making the net kinetic mixing term equal to zero. This mechanism is analogous to the 
anomaly cancellation one where the GS mechanism can be generated by an anomaly free theory with 
some heavy chiral fermion integrated out of the mass spectrum [24] . 

7 We stress that the ips — ^A sector presents a different parameters choice with respect to [39]- [41], 
where we realized a scenario in which the mixing between tjjs and was suppressed. 
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just considering the superpotential (3), the gaugino masses (19) and by reminding that 
$r gets a vev i>_ different from zero, while v + = 0. 

Finally there are also corrections coming from the anomalous axino couplings: F-term 
couplings of the type b aa (F s ) A a A a and D-term couplings of the type b aa ^s^a(D a ), and cor- 
rections coming from the superpotential term e~ kS H u H d + h.c. However such corrections 
are always subdominant and thus we neglect them with very good approximation. 

We assume the lightest supersymmetric particle (LSP) in our model comes from the 
neutralino sector. In Sec. 8 we show the parameter regions in which this holds true. In 
order to ensure that the neutralino is the LSP we keep fixed the gravitino mass m 3 / 2 ~ 
O(TeV) in the limit M P -> oo. 



7 Sfermion masses 

The sfermion masses receive several contributions. As we have seen in Sec. 3 the leading 
contribution comes from the induced soft masses (17). But there are further contribu- 
tions. We have MSSM-like contributions: F-term corrections proportional to the Yukawa 
couplings, Dy and D 2 term correction from the Higgs sector. Moreover there are Da term 
corrections from the Higgs and _ sector. As an aside, the appearance of such terms in 
the low-energy action, given our assumption My A > 1 TeV, can be understood in terms 
of quantum corrections to Kahler potential [34]. Considering the first two families we 
neglect the corresponding Yukawa couplings (the so called third family approximation). 
In this approximation the mass eigenvalues are given by 

2 2 2 1 / 1 2 21 

///,.. ~ 111.. =m; ' ■ ■ 1 •■ ■ • 



U L — " U C L - Q 



{l g2y ~ 92 ) + qQrh2 ° A (74) 

m L - m L = m h-9 Y ^- + Qu^m 2 DA (75) 

f\ \ Aw 2 

m \ L ~ m L = m % + ( + 9 2 2 J -y + qgrri 2 DA (76) 

A?; 2 

2 2 2,2^^1 ~2 

m d R ~ m ~SB = m Dc + 9Y-^ + QD^m DA (77) 

A?; 2 

2 2 2 / 2 1 2 \ 1 ~ 2 r>7o\ 

m u e = m v, = m L-{9 Y + 9 2 )-^ + ( lLm DA (78) 

Av 2 

m l L ~ m l L = m L ~ (9y - 9l) -y + qLm 2 DA (79) 

Av 2 

m L - m L = m % + + <lEcm 2 D A (80) 

The first terms on the right hand side m% ~ ~ ~ ~ are the corresponding soft masses 
(17), the second terms are the D Y ,2 contributions with At> 2 = v 2 u — v 2 d = — f 2 c 2/ 3, while the 
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last terms are the Da corrections given by 

1 



™d a = o {<te« v l + q Hd v 2 d - vi) (81) 



There is an approximated degeneracy between the sfermions with the same charges. 
The mass matrix for the third family sfermions is parametrized as 

^ (m''m''\ 

where the off-diagonal terms are generated by F-term corrections proportional to the 
Yukawa couplings. The stop mass matrix elements are 

- 2 / 1 \ A?; 2 

Jlrft _ ™2 i _2 / 1 „2 „2 \ L ^ u i „ ^,2 



ll = m t+ m Q+ [ g^y - 92 J -y + <?q™d a 

6~ 



m rr = m t+ m lc - 9y^tt- + Qu°rh 2 DA 



The sbottom mass matrix elements are 



Ml R 2 = -fim^ 1 (83) 



rb 2 _ _2 , _2 , f^„2 , „2^ ^ 2 



= ™& + ™J + (^y + 92 j -y + <?Q™D A 
M| R 2 = -//m 6 ^ (84) 



The stau mass matrix elements are 

2 



At; 2 



Mr 2 2 , 2 / 2 2\ ^ u , ~ 2 

ll =m T + m L - {g Y - g 2 ) — + q L m DA 

» *-f 2 2 o 2 ^ 2 ~ 2 

= -fJ.m T tf} (85) 

The tau sneutrino mass is 

Aw 2 

m i =m \~ (9y + #2) -g— + QLm 2 DA (86) 

where m t , m;, and m r are the masses of the corresponding standard fermions (i.e. further 
F-term contributions proportional to the Yukawa couplings). The structure of the diagonal 
terms of (82) is the same as in eq. (74)-(80): soft masses, MSSM D-term contribution and 
Da term correction. Furthermore we stress that there is a mass degeneracy between the 
three sneutrinos £ e ,/i,r since the soft masses (17) are flavor blind. 
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Figure 1: Higgs couplings bounds. The yellow spot represents our charge choice. 

8 Phenomenology 

In the following we derive the phenomenological consequences of our scenario. Following 
our assumption of having a mass parameter for the anomalous U(1)a just slightly above 
the TeV scale, we fix My A = 10 TeV. The mass scale in the gaugino sector A is set to be 
0(M Va ). 

8.1 Charge Bounds 

The model parameter space can in principle be constrained by precision EW measure- 
ments [35]. However, since My A = 10 TeV every value of gAQH u and gAQH d is allowed by 
EW precision data if \gAQH u \, \9AQH d \ J$ 0.1. So the only relevant constraints are (27) and 
(54), that are plotted respectively with a red and a blue region, in Fig. 1 in the plane 
{QAQHu, 9AqH d )- 

8.2 Free Parameters 

Here we discuss which parameters remain free in our model after all the constraints 
discussed in the previous sections are imposed. Our choice for the Higgs U(1)a charges 
corresponds to the yellow spot in Fig. 1 

g A = 0.1 M Va = 10 TeV 

q Hd = -(1/3) q Hu = -(2/5) (87) 
In order to fix the remaining parameters ((a), (Fs), A, m, Qb) we assume v ~ 246 
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GeV and then we choose some benchmark value for g B and V- in the U(1)b sector 8 : 

A) g B = 0.4 v. = 5 TeV (88) 

B) g B = 0.1 w_=4TeV (89) 

The next step is to solve the minima conditions (41)-(43) determining (Fs), A, m as 
function of (a). In the limit in which v 2 <C My A (a),v 2 _, we get 



1 2 (">9A('?H d +'3g tl ) 

A 2 ~ 4e 



fi'A (^9a v - - 2(a)M VA ) (sec{2(3){q Hd - q H J + q Hd + q H J 



<*s> ^ -e ^ 77- ^— 2(a)M^-^ 2 



m 



- ^(«)M yA -l[(^ + ^)^] (90) 

In Appendix B we report the exact formulae. Thus the only remaining free parameters 
are tp and (a) and we perform the following analysis of the mass spectrum as a function 
of tp and (a). A lower bound on (a) as a function of i>_ can be obtained, given the 
approximation (45), from the equation (48) 

, , H 2 + (i/2) (A + <A)vi .... 

<C,) = 9^ (91) 

where we used the relation 

ml- = -ml = -(a)g A M VA (92) 
Thus the lower bound on (a) is obtained simply by setting \m\ =0, 

{ >m - 3aM Va [MI 
The condition (a) > (a) m must hold since otherwise we would have a massless scalar 
field in the spectrum (see eq. (55)). Another lower bound, (a)&, can be obtained from 
the condition (39), by solving the minima conditions (41)- (43) and by substituting the 
corresponding (F s ), A and m values (90). The resulting lower bound can be expressed as 

(a) > max[(a) m , (a) b ] (94) 

No upper bound can be imposed, hence we decide to perform our analysis by considering 
(a) < 100 TeV. 

The parameters A and (Fs) are of a particular phenomenological importance since 
they appear in the /x and b terms (see eqs. (32) and (33)). In the case A, \i is in the range 
(900,6000) GeV and Vb is in the range (50, 1200) GeV while in the case B, /i is in the 
range (500, 6000) GeV and Vb is in the range (25, 1200) GeV. These values are in the 
right range to solve the /i-problem. 



3 We remind that v + — (see Section (4)). 
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8.3 Mass spectrum 



A) With such choice the gauge vector sector is completely fixed up to a tp dependence. 
Anyway even such a dependence can be safely ignored with a very good approximation 
in the new gauge sector since the mixing is strongly suppressed. So for each tp value we 
have 

M Zl ~ 10 TeV (95) 
M Z2 ~ 2 TeV (96) 

where with Z 1 we denote the V^-like vector . 



B) As in the previous case, we just give the Z 12 masses 

M Zl ^ 10 TeV (97) 
M Z2 ~ 400 GeV (98) 

where as in the previous case Z\ is Vyi-like. 

We will not give the exact values of the Zq mass. It is enough for our purposes to 
know that they are compatible with the bounds of Section 8.1. Both case A and B are 
compatible with CDF bounds about Z 1 direct production [42]. 

Recent LHC data have restricted the most probable range for the Higgs particle mass 
to be [115.5,131] GeV (ATLAS) [47] and [114.5,127] (CMS) [46]. Moreover, there are 
hints observed by both CMS and ATLAS of an excess of events that might correspond 
to decays of a Higgs particle with a mass in a range close to 125 GeV. So, in Fig. 2 
and 3 we give region plots showing the allowed values of (a) and tp for case A (B) and 
A/^/c = (5)10 TeV. The red region is the one in which M^,^ loop e [124, 126] GeV where 
the h° mass is computed considering 1-loop corrections. Since it turns out that the top 
squarks have small mixing angle and considering the limit M^o ^> Mz , we have [36] 

M ^li-ioop - M h°Loc + ^ s iyt m t ln {m.m-JmVj 

3 TTt^ 

- M ^lt r ec+^^ ln K^ 2 /^ 2 ) (") 



where Mho\ tree is the tree-level h° mass and we used m t = y t v u /2 = y t vsp/2. There is 
an approximated inverse correlation between (a) and tp in the ho mass allowed region 
because the 1-loop correction in (99) increases for increasing values of (a) or tp. The ho 
mass allowed region is almost the same for case A and B because of two reasons 



i. the mixing with ± is suppressed 
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CASE A, A c =5, m->30& GeV CASE A, A c =10, m s ->308 GeV 




20 40 60 80 100 20 40 60 80 100 

<a> (TeV) <ff> (TeV) 



CASE B, A c =5, m ->308 GeV CASE B, A,.=10, m->308 GeV 




20 40 60 80 100 20 40 60 80 100 

<a> (TeV) <ff> (TeV) 



Figure 2: Allowed (a) and tan/3 values for case A (up) and case B (down), A c = 5 (left) 
and A c = 10 (right). The red region is the one in which -^o| lloo £ [124,126] GeV, 
the magenta region is the one in which M^o^ loop G [114.5, 131] GeV and the blue region 
satisfies all the mass bounds on the sparticles (from PDG) and requires a neutralino LSP. 
The yellow dots are our benchmark points. 

ii. the parameters rhh u ,rhh d , fi,b in the scalar potential (52) are ruled by the square 
mass parameters gA{ot)My A and (gAV-) 2 and the first one turns out to be dominant. 

The magenta region satisfies a milder constraint on the light Higgs boson: M^ 1 1 loop E 
[114.5, 131] GeV. In order to be more conservative we imposed the joint constraints of 
ATLAS and CMS. 

The blue region satisfies all the mass bounds on the sparticles and requires a neu- 
tralino LSP. We considered two possibilities: one more optimistic (Fig. 2) using the PDG 
bounds [43, 44] and one more conservative (Fig. 3) using recent LHC data [45]. The 
combination of the gluino mass bound with a neutralino LSP is a strong constraint that 
reduces drastically the allowed parameter space. In some cases there is not even a blue 
region, which means that we cannot satisfy simultaneously all the mass bounds and have 
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CASE A, A ( .=5,m ->1.1 TeV 




CASE A, A c =10, m->l.l TeV 



20 40 60 80 

<ff> (TeV) 
CASE B, A =5, m , ; >1.1 TeV 





100 20 40 60 80 

<ff> (TeV) 
CASE B, A,. = 10, m , ; >1.1 TeV 




<«>(TeV) 



<«> (TeV) 



Figure 3: Allowed (a) and tan/3 values for case A (up) and case B (down), A c = 5 (left) 
and A c = 10 (right). The red region is the one in which -^o| lloo £ [124,126] GeV, 
the magenta region is the one in which M^o^ loop G [114.5, 131] GeV and the blue region 
satisfies all the mass bounds on the sparticles (from preliminary LHC data) and requires 
a neutralino LSP. The yellow dots are our benchmark points. 



a neutralino LSP, so they are completely ruled out. When the gluino mass bound is 
from PDG, Case A is allowed, otherwise it is completely ruled out, and only case B for 
A/y/c = 10 TeV presents allowed regions. We notice that case A favors low (a) values, 
while case B favors big (a) values. For every allowed case we choose a benchmark point 
(yellow spots in Fig. 2 and Fig. 3) 



i. case A, A c = 5, (a) = 3 TeV and tp = 50 so that Mho\ v 



loop 



121.6 GeV 



ii. case A, A c = 10, (a) = 5 TeV and tp = 10 so that M h o\ 1Aoop ~ 124.7 GeV. 

iii. case B, A c = 10, (a) = 4.5 TeV and tp = 50 so that M h o| l loop ~ 125.1 GeV. 

-loop 



iv. case B, A c = 10, (a) = 50 TeV and tp = 2.5 so that M h o\ 1Jn _ ~ 130.1 GeV. 
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CASE A, f j8 =50, A c =5, <a>=3 TeV, Mass spectrum 
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CASE A, ^=10, A t .=10, <ff>=5 TeV, Mass spectrum 
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Figure 4: Mass spectrum, case A, A c = 5, (a) = 0.3 TeV and tp = 50 (left), A c = 10, 
(a) = 0.5 TeV and tp = 10 (right). 



CASE B, ^=50, A c =10, <o->=4.5 TeV, Mass spectrum 
3.0 [ 



2.5 - 



CASE B, r«=2.5, A c =10, <a>=50 TeV, Mass spectrum 
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Figure 5: Mass spectrum, case B, A c = 10, (a) = 0.45 TeV and tp = 50 (left), (a) 
TeV and tp = 2.5 (right). 
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and we give the full mass spectrum in Fig. 4 and in Fig. 5. 
All the benchmark points share some common features 

• the LSP is the lightest neutralino of the new sector: in case A it is a combination 
of ± and \ B while in case B is almost a pure A#- 

• an approximated mass degeneracy of H°, A and H ± holds, and their masses satisfy 
the bounds of [38, 48]. 

• the lightest sleptons is a sneutrino, except for tp = 50 when it is f\ 

• the lightest squark is Ul, except for tp = 50 when it is b\ 

• the first and second family left-handed squarks / sleptons are likely to be lighter than 
their right-handed counterparts. This is at odds with the usual MSSM cases [36]. 

• Ci(2) is close in mass with N^f? M . C 2 and N^ SSM are heavier than all sfermions. 

• the gluino is close in mass to C\ and A^^ SSM which are gaugino-like. Moreover 
it is lighter than all the squarks except for point i). So it turns out to be long 
lived, specially in case B where the approximated mass degeneracy involves also the 
LSP. Long lived gluinos bind with SM quarks and gluons from the vacuum during 
the hadronisation process, and produce R-hadrons. R-hadrons are among the most 
interesting searches at LHC. Anyway we will come back to this point with a more 
detailed study in a forthcoming paper. 

• there is an approximated mass degeneracy between e R and u R because using the 
charge constraints (25) and (87) we get g^c = 3 and que ~ 2.9. 

• m^- < except for point i) 

Case B points deserve some more comments. 

(j) + and iV£| w are out of the plot of point iv) because they are heavier than 6 TeV. Z 2 
is among the lightest not SM particle, so it can decay only into SM particles, because of 
energy and R-parity conservation. So Z 2 is long lived, because SM particles are coupled 
to Z 2 only through the suppressed Va,b mixing or through the Higgs scalars which present 
a tiny mixing with ± . 

It is not an easy task to compare the resulting spectrum we get for our model with 
those related to the rich zoology of supersymmetry breaking scenarios. It is worth to 
stress anyway that the two representative spectrums showed in Fig. 5 which encode the 
key features of our scenarios listed above are not reproduced in any of the benchmark 
points showed in [49]. 
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9 Conclusions 



In this paper we presented a viable mechanism to generate soft supersymmetry breaking 
terms in the framework of a minimal supersymmetric anomalous extension of the SM. 
The crucial ingredient is a non perturbative term in the superpotential (3) which couples 
the Stiickelberg field S to the Higgs sector. This term is related to the generation of a 
suitable \x and b terms (see Eq. (32) and (33)) in the low energy effective action when the 
Stiickelberg gets vev. We argued about the origin of this term from an exotic instanton in 
an intersecting D-brane setup. We computed the spectrum of our model as a function of 
the saxion vev (a) and for different choices of the remaining free parameters. We checked 
our results against known phenomenological bounds, namely current lower bounds on 
the mass of the scalar and fermionic superpartners. We analyzed a scenario in which 
the anomalous sector is the source of the soft supersymmetry breaking terms while the 
corresponding vector and Stiickelberg multiplets are not present in the low energy effective 
action. For what concern the non anomalous sector we took into account two different 
cases (dubbed case A and case B). 

As we stated in Sec. 8, by applying some phenomenological constraints we were able 
to derive some bounds on the saxion vev (a), which is the relevant parameter setting 
the mass scale of the scalars. The strongest constraints on (a) and tp comes from the 
combined requirement of -^o| lloo £ [124,126] GeV or ([114.5,131] GeV), a neutralino 
LSP and that all mass bounds (specially the gluino one) are fulfilled. In Fig. 2 (pre- 
LHC bounds) and 3 (preliminary LHC bounds) we summarize the allowed regions for 
(a). In the first case, by requiring a phenomenological appealing neutralino LSP, we get 
an allowed (a) of few TeV up to 10 TeV for the A and B scenarios respectively. In the 
second case (preliminary LHC bounds) we get that only the B scenario is allowed with 
(a) > 5 TeV. These results can be seen as a bound that a concrete D-brane model has to 
satisfy. We deserve this analysis for future work. 

In Fig. 5 we explicitly showed two benchmark mass spectrums for our model with (a) 
and tp which fulfill the above bounds. The cases shared different peculiar features: the 
LSP is the lightest neutralino of the new sector, there is a near mass degeneracy between 
H°, A and H ± , and between e R and u R , the lightest sleptons is a sneutrino except for 
tp = 50 when it is stau, the lightest squark is a u L except for tp = 50 when it is a 
sbottom, the first and second family left-handed squarks/sleptons are typically lighter 
than their right-handed counterparts. Moreover in case B the gluino is long lived and 
can produce R-hadrons. It turns out that these features are not reproduced in any of the 
widely studied benchmark points presented in [49]. 
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A Anomalous Lagrangians 



The Lagrangian involved in the anomaly cancellation procedure is 



1 



(S + S* + 2M VA V A y 



e 2 e 2 



(100) 



9 2 a b aa STi {W a W a ) + g Y g A b ¥A SW Y W A 



+ h.c. 
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where the index a = A,B,Y,2,3 runs over the U(1) A , U(1) B , U(1) Y , SU(2) and SU(3) 
gauge groups respectively, and the constants b ab are fixed by the anomaly cancellation. 

Since we have only one anomalous U (1) we can avoid the use of GCS terms, distributing 
the anomalies only on the U(1) A vertices. So we have 

9aAaa yy 9aAyy 



b AA = 
b 33 = 



96it 2 M Va 
9aA 33 



b YY = 
b YA = 



b 22 = 



32tt 2 M Va 
9aAya 
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9aA 
16% 2 M Va 



1Qtt 2 M Va 32it 2 M Va 
where the ^4's are the corresponding anomalies 

Aaa = -l0qH d 3 -9qH d 2 (qL + 3q Q )-9q Hd (q L 2 + 3q Q 2 ) 
-7q Hu 3 - 27q H 2 q Q - 27q Hu q Q 2 + 3q L 3 

A YY = --(7q Hd + 7q Hu + 3q L + 9q Q ) 

A22 = ^(qH d + qn u + 3q L + 9q Q ) 

A33 = -^(qH d + qH u ) 

A Y a = ^qn d 2 + QqH d (qL + qq) ~ qH u (5q Hu + l2q Q ) 



(101) 



(102) 
(103) 

(104) 

(105) 
(106) 
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where we used the constraints (22). Imposing the conditions (25) we get 

Aaa = ^(-U68q Hd 3 + l776q Hd 2 q Hu - 996q Hd q Hu 2 + 53q H J) (107) 

Ayy = -^(qH d + qH u ) (108) 

A22 = -\(QH d + QH u ) (109) 

^33 = -\{qH d + q Hu ) (110) 

A YA = (111) 

We remind that (111) is not a consequence of (25), but rather (25) is a consequence of 
imposing (111) in order to cancel the U(l)y — U(1) A kinetic mixing. 



B Exact fixed parameters 



In this Appendix we give the exact values for the (i*s), A, m parameters determined in 
section 8.2. Solving the minima conditions (41)-(43), we get 



e 2{a)g A (q Hd +q Hu )/M VA 
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x 



(112) 



x 



- g A sec(2(3)(q Hd - q H J (8(a)M VA + g A v 2 (q Hd + q Hu ) (cos(4/3) + 3) - Ag A v 2 _ ) 



-8(a}g A M VA (q Hd + q Hu ) - 2v 2 (2g 2 A {q H 2 + q H 2 ) + g\ + g 2 2 ) + 4g 2 A v 2 _(q Hd + q Hu 



,( a )9A( 



(iH d + q H u )/M VA x M FA tan(2^) 



x 



8gA(qH d + qn u )\ 
9A{qu d ~ qn u ) (4(a) M Va + g A v 2 (q Hd + qH, 

cos(2/3) (g 2 A (q Hd - q Hu ) 2 + g\ + g 2 2 ) 



x 



2g A v 2 _) 



H 2 = g A {ot)M VA - - [(g 2 A + g 2 B ) v 2 _ + g\v 2 (q Hd c} + q Hu s})] 
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